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ABSTRACT 

A study of the semi-analytical approach to four-fermion production in e~^e~ annihila- 
tion is presented. We classify all possible four-fermion final states and present results of 
new calculations for the 'basic' processes with the WW, ZZ, and ZH off-shell production 
together with some examples of 'background' processes. The Initial State Radiative cor- 
rections are included for the basic processes. Several numerical examples are given in the 
energy range from LEP 2 up to ^/s = 1 TeV. 
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1 INTRODUCTION 



While at LEP 1 the basic process is two-fermion production via a single Z-boson (or 
photon) exchange, at LEP 2 the typical process will be the 'double-resonance' production 
of four-fermion final states. LEP 2 will operate just in the threshold region of W- or Z-pair 
production and, if the Higgs-boson mass fulfils 60 < Mh < 100 GeV, of ZH production. 
Already in the tree level approximation, the double-resonance physics is much richer and 
much more interesting and complex than the single Z-boson production. 

A specific four-fermion final state can be produced by many Feynman diagrams with 
many possible virtual states, including all the carriers of fermion interactions in the Stan- 
dard Model: 'y,Z,W'^, g, and H. We will distinguish between 'basic' diagrams, which 
contain two potentially resonating virtual states {W, Z, H) in the s-channel and 'back- 
ground' diagrams, which are just the rest. The general topology of the basic diagrams 
is shown in figure 1. The contribution of background diagrams to a given final state 
is usually suppressed. Some background diagrams for IV-pair production are given in 
figure 2. 

The analytical result for the Born on-shell VT-pair production has been known for long 
in the literature The off-shell case was treated in [^. A calculation of off-shell Z-pair 
production was done in 0]. 

In this paper we present semi-analytical results for off-shell production of bosonic pairs 
(WW, ZZ, ZH) including universal lowest-order Initial State Radiative (ISR) corrections 
with soft-photon exponentiation. For the case of W^-pairs, we also present results with 
complete 0{a) ISR corrections and some examples of four-fermion background processes. 

The paper is organized as follows: In the next section, we classify the four-fermion 
production processes in e~^e~ annihilation. In section 3 we introduce the notations and 
present the formulae for the basic processes. In section 4, the contributions of the back- 
ground diagrams to the simplest final state are characterized. The ISR corrections are 
described in section 5. In section 6, we present concluding remarks and prospects. 

2 A CLASSIFICATION OF FOUR-FERMION 
PROCESSES 

In this section we will classify the four-fermion production in the Standard Model i. The 
number of Feynman diagrams depends crucially on the final state. In general all possible 
final states can be subdivided into two classes. 

The first class comprises production of up (anti-up) and anti-down (down) fermion 
pairs, 

[U, A) + [D, Uj) , 

where i, j are generation indices. The final states produced via virtual H^-pairs belong to 
this class. We will call these 'CC type final states. The second class is the production of 

"'^Some approximate results on the analytical treatment of off-shell VF-pair production near threshold 
can be found in . 

The classification is done with the use of CompHEP (5|. 
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Figure 1: The basic contributions to off-shell boson-pair production: crayfish and crab; 
B = W^,Z,H. 

two fermion-antifermion pairs, 

{fi fi) + (/,• /,) J-U, D. 

As it is produced via a pair of two virtual neutral vector bosons we will call this a final 
state of the 'NC type. Obviously these two classes overlap for certain final states. 

The number of Feynman diagrams in the 'CC class is shown in table 1. Three different 
cases occur in the table: 
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Table 1: Number of Feynman diagrams contributing to the production of 'CC type final 
states. 



(i) The two produced fermion pairs are different {i j) and the final state does not 
contain an (numbers in boldface). For this case, the number of diagrams varies 
between 9 and 11, depending on the final state's neutrino content. The background 
diagrams for this simplest case are shown in figure 2. 

(ii) The four- fermion final state contains one ePg- or ei/g-pair (roman numbers); the 
number of diagrams grows to 18, 19 or 20, due to the additional t-channel exchange 
diagrams. 

(iii) Two mutually charge-conjugated fermion pairs (i = j) are produced {italic num- 
bers). Here, the diagrams may contain neutral-boson (Z, 7, H, gluon) exchanges. 
One should emphasize that this overlaps the 'NC classification. 

For the final states corresponding to the 'NC class the number of Feynman diagrams 
is presented in table 2: 
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Figure 2: Background contributions to off-shell VT-pair production: up and down 
reindeers. 

(i) The simplest case (numbers in boldface) does not contain electrons or identical 
fermionsQ. 

(ii) With identical fermions / (/ 7^ e, z/e), the number of diagrams (in typewriter) 
grows drastically, since it is necessary to satisfy the Pauli principle (i.e. to antisym- 
metrize the amplitude). 

(iii) The numbers in romans correspond to the final states that include / = e, z/g except 
those covered by item (iv). The large number of diagrams here is due to additional 
t-channel diagrams. 

(iv) The numbers in italics correspond to final states that are also present in table 1, 
case (iii). The basic diagrams proceed via both WW- and ZZ -exchanges. 

So far we have investigated semi-analytically only the two simplest cases of the 'CC 
and the 'NC classifications. In section 4 we will present the production of two different 
fermion pairs: 

e+e- ^ {W+W-,W^ff') f'Mh, 

fijj 7^ e, Z/e . 

An example of the 'NC type process with 24 Feynman diagrams (case (i) of table 2) will be 
discussed in another contribution to this conference p. For the more involved processes, 
especially for case (ii) of table 1 with final states containing an eue pair, investigations 
have been started. We expect that the topologies indicated with roman numbers may 
also be treated by our method. 

Another classification of four-fermion processes was presented recently . This refer- 
ence separates leptonic, semileptonic, and hadronic four-fermion final states. For leptonic 
processes they agree with our number of Feynman diagrams. For semileptonic processes 
agreement is found in all cases but one, namely qqVf.Vf,. The authors of [0] have 19 dia- 
grams instead of our 21 for this case. The explanation is that for this process a H^-fusion 
diagram exists, where two ly-bosons produce a qq state with a virtual q' in the t-channel. 

^We exclude the Higgs-boson exchange diagrams from the classification in the tables. 
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Table 2: Number of Feynman diagrams contributing to the 'NC type production of two 
fermions pairs. 

For a given g, say d, q' may be u,c or t (due to Cabibbo-Kobayashi-Maskawa mixing). 
In our classification we count all three diagrams, while in this is counted only once. 
Finally, for hadronic processes we agree only in one of seven cases, namely in our example 
with 11 diagrams. The reasons for these differences are twofold. Firstly, once more, quark 
mixing is neglected in 0. Secondly, we count gluon-exchange diagrams while this is not 
done in ||^. Taking these differences into account, both classifications agree. 

At the end of this section we introduce some notations. Background diagrams may 
contain one resonating virtual s-channel state or none (non- resonating background). We 
will denote by a sub-index n the total number of resonating s-channel propagators in a 
separate contribution cr„ to the cross-section of four-fermion production, 

e+e" ^ {W-^W-,ZZ,ZH,...) ^4f. 

Therefore, the basic contributions carry sub-index 4, the background-basic interferences 
3 or 2, and the pure background contributions have sub-index 2, 1 or 0. With appropriate 
kinematical cuts, certain resonating states may be selected and contributions with indices 
4, 3, 2, 1, and should be hierarchically less and less important. 



3 BASIC CROSS-SECTIONS 

Here we present the basic cross-section formulae for the three cases described by the 
diagrams of figure 1 with WW, ZZ, and ZH intermediate states. They all may be 
expressed by twofold convolutions of a hard-scattering off-shell cross-section with Breit- 
Wigner density functions. For the WW case, this representation was invented in 0: 



(v^-v^) 
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cr^^(s) = J dsipwisi) J ds2Pw{.S2) 


X af^(s;si,S2), (1) 
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where 



Pw[Si 



SiTwisi) X BR(i) 



TT \si - + iy/TiVwisi 



(2) 



is the Breit-Wigner density function originating from the s-channel propagators 
and BR(z) is the corresponding branching ratio. Similar densities pz and pn associ- 
ated with Z {H) s-channel exchanges can be obtained by the replacements Mw, 
Mz, Tz] {Mh, Th)- They are normalized so that 



Pb{s) 
Further, 



6{s - Mb) X BR{i), B = W,Z, H. 



DTTV 2 J 



(3) 



(4) 



is the off-shell (s-dependent) VT-width and Nc{f)=l{3) for leptons(quarks). The sum 
in (^) extends over all open fermion channels. At extremely high energies, additionally 
opening channels may substantially contribute to rw{s). 

The cross-section {s, si, S2) contains six pieces but is described by only three 
functions ^4(5; Si, S2), which are different for the s- and t-channel and the st-interference: 



WW, 



S] Si, S2) 
2 \2 



■ins 



(c^^ + c^z + czz)QI{s] Si, S2) 
+ {cu'y + c^z)Gf{s;si,S2) 

+ C^uQlis] Si, S2) 



(5) 



The coefficients Ca/3 consist of Z(7)ee- and Z(7)l^l4^-couplings and 7, Z-propagator ratios: 
c^-y = 8s^Q1, 



c^z = 4:s'^Ve\Qe\'^e 



s-Ml + zMzTz{s) 



czz 

CuZ 



S-MI + iMzTz{s) 



-4s^ \Qe 



s-Ml + iMzTz{s) 



with Qe = ~ 
functions are 



(6) 



l,ae = l,fe = 1 — 4s^,s^ = sin 9^. The three irreducible kinematical 



Ql{s]Si,S2) 



S^SiS2 



- + 2s(si + S2) + 2siS2 
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gf{s;si,S2] 



(X 



s'^SiS2 \ 6 



s + ll(si + S2)] 



+ 2s{sj + 3siS2 + si) - 2{sl + si 

- 4siS2 [s{Si + S2) + S1S2] A 



^4(s;si,S2) 



AV2 



- + 2s S1 + S2 



with 



and 



SS1S2 

- 8S1S2 + 4SiS2(s - Si - S2)A 

= A(s;si,S2) 

= s'^ + sl + si — 2ssi — 2siS2 — 2S2S 



^ , . 1s-si-S2 + \/A 

A(s;si,S2) = ^ In 7= 

vA s — si — S2 — vA 



(7) 



(8) 



(9) 



In the hmit @, the on-shell W^-pair production cross-section is obtained as given in 
The corresponding set of formulae for basic off-shell ZZ production is: 



a^^{s\ 



dsipz{si) / ds2Pz{s2) 




ZZf 



X at^{s;si,S2). (10) 
The off-shell Z-boson width, which enters the definition of pz, has the form 
G„ Ml 

Vs. 
' / 



Tzis] 



247r 



(11) 



The cross-section erf ■^(s, Si, S2) is extremely compact and can be described by only one 
function Q^^'^{s] Si, S2), which is the sum of three others: 



Qi{s]Si,S2) = ^4(3; Si, S2) + ^4(5; Si, S2) 

+ GTis;si,S2). 



(12) 



Here the functions ^4(s;si,S2) = ^4(s;s2,Si) correspond to u- and t-channel diagrams 
and ^4"(s, Si, S2) describes the tw-interference. The cross-section erf ^ is given by: 



a^^{s; si, S2) = 
^^l^{vt + Qv!al + at)gt'-{s,s„S2) 



(13) 
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with the kinematical function 



't+U| 



S', Si, S2) 



+ {Sl + 32^ 
S - Si - S2 



Finally, the basic off-shell ZH cross-section is: 

dsipu{si) \ ds2pz{s2] 



X af"{s]Si,S2) 



Below the threshold of the decay H —>■ W~^W , the off-shell width has the form 



TMis) 



G 



(14) 



(15) 



The cross-section af ^(s; si, S2) is given by 



7 T-f / \ 

0-4 (s;si,S2) 



X 



967rs s 
2 



Q^Hs] Sl, 82). 



S-MI + iMzTz{s) 
Again, it contains only one kinematical function, namely Q^\s] Si, S2): 



(16) 



Q^\s] 81,82) 



AV2 

-^{\ + 12882) 



(17) 



The energy dependences of the off-shell WW , ZZ, and ZH basic cross-sections are 
presented in figures 3-5. For comparison we also show on-shell cross-sections for all 
three cases. With respect to the on-shell case, off-shell cross-sections are substantially 
reduced in the threshold region and develop tails at high energies so that boson widths 
cannot be neglected for s ^ {Mb-^ + M^J^. Going off-shell, the cross-section peaks are 
shifted to higher energies. Using constant widths, i.e. y/sTsis) — > MbTb, corresponds 
to redefinitions of the boson masses ||: Mb ^ Mb = Mb + \T%/Mb. For the W this 
results in Mw ~ Mw + 26 MeV [Q. That lowers the cross-section around threshold by 
at most ~ —1.7%, while for -y/s > 180 GeV the effect is small and positive (< 0.1%). 
Numerical results were obtained with the Fortran program GENTLE |TD . 



4 BACKGROUND 

In this section we briefly sketch our results for the semi-analytical treatment of the pro- 
cess (H) , the production of four different fermions excluding electron or electron neutrino 
(case (i) in table 1). This is described by the basic diagrams of figure 1 and by the 
eight background diagrams of figure 2. In the unitary gauge, there are no other dia- 



grams for this particular final state. All results were obtained with the help of FORM |]Tl 
SCHOONSCHIP ra, and CompHEP H. 
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Figure 3: The basic H^-pair production. 



Figure 4: The basic Z-pair production. 
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Figure 5: The basic ZH production. 



After the analytical integration over five angular variables, one arrives at a doubly- 
convoluted representation for the cross-section: 



are interference contributions. 
We obtained explicit representations for the cross-sections crt^ {s; Si, S2) and <72^ {s; Si, S2) 
in terms of seven new kinematical functions (three functions and four functions ^f), 
coupling constants and propagator ratios (similar to (^). The formulae for the cross- 
sections themselves are rather appealing, while only three of the kinematical functions 
are of the same compactness as those of the basic processes. Four interference functions 
(two ^3 and two ) could be written only with the aid of a cumbersome polynomial 
presentation of the following type: 




1 



^3 (s. Si, S2) = ^/XYI [■5iA(s; si, S2)Y 



i,j=0 



X [5152^4(5; si, S2)y VIAs] Si, S2), 



(19) 



1 



^2 (s, -51, S2) = X! [■5iA(s; si, 82)]' 



i,j=0 



X [s2C3{s;s2,Si)yV^As; 81,82), 



(20) 
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3 /'„ „ \n-l 



E 

n=l 

with £4(5; Si, S2) given by (H) and 



+ E ^^^1^— (21) 



si, S2) = ^ In — Ty- (22) 

VA s + Si - S2 - V A 

In (ig) the a stands for f/, D, and in for UD, UU. The p^(u), = 0, 1, 2, 3 in (grj) 
are polynomials of order n in s, si,S2. A complete analytical result will be presented 
elsewhere \T3\. Here we restrict ourselves to a comment and some numerical results. As 



may be seen from (jT9|)-(|2TD, the cumbersome interference functions contain inverse powers 
of A (up to the third power), which vanish at the upper limit of integration over S2- This 
is a typical example of so-called kinematical singularities. Expanding £3^4 in Taylor series 
in A, one may see that all these inverse powers cancel and the cross-section has a proper 
threshold behaviour. However, these kinematical singularities may create complications 
for numerical calculations. 

In figure 6, we present the ratios (basic-|-background) /(basic) , 

for three different channels with 11, 10, and 9 diagrams respectively (see case (i) of table 
1) as functions of a/s. As is seen from the figure, the background contributions for these 
processes are relatively small, especially at LEP 2 energies. Below threshold of the W- 
production and at high energies, the relative contribution of such background increases 
and reaches ~ 2%. 



5 INITIAL STATE RADIATIVE CORRECTIONS 
IN FOUR-FERMION PRODUCTION PROCESSES 

Since the background contribution is comparatively small, it is quite reasonable to re- 
strict oneself to the basic diagrams when calculating the complete lowest-order ISR QED 
corrections. 

Let us recall that the ISR corrections are known to be dominating in single Z-resonance 
production. A similar property holds for four-fermion processes at LEP 2, which will 
proceed near the corresponding i?iS2-thresholds and the intermediate-state bosons will 
be nearly at rest fl^ . 



The ISR corrections for a cross-section described by s-channel diagrams with 7- and 
Z-exchanges (i.e. including the background terms) may be presented by a universal 
formula 



2 
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Figure 6: The (basic+background) /basic ratio R for off-shell Vl^-pair production with 
k ^ e. 



ds' r 



[3,v^^-\l + S) + H 



(24) 



where v = 1 — s' /s. The soft plus virtual photon part 5* and the hard part H{s'/s) are 
given by 



S 

H{s'/s) 




+ + 0{a') 



(25) 
(26) 



and = 2a/iT[ln{s/ml) — 1]. 

Equation (|2^) may be directly applied to the case -B1-B2 = ZH and to the s-channel 
contribution of -B1-B2 = WW. We have rederived by explicit calculations that ( |2^ may 
be obtained straightforwardly from the usual vertex and bremsstrahlung QED Feynman 
diagrams after up to seven sequential angular integrations (five for the vertex part and 
seven for bremsstrahlung). The situation becomes much more complicated if t- and u- 
channel exchange diagrams are involved. Here we face two kinds of problems. The 
factorized form (^4]) is no longer valid for the squared t- and w-channel diagrams and 
for st-,su-, and t-u-interferences. One of the reasons is the angular dependence of t- and 
M-channel propagators. This leads to the appearance of additional non-factorizable (non- 
universal) QED corrections. 
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For the ZZ process this is the only problem, although technical complications arise 
due to additional Feynman diagrams with real and virtual photons attached to the virtual 
electron line in t- and u-channels. For the WW basic diagrams another problem persists. 
There, the electric charge flows from the initial state electron through an intermediate 
VT-boson to a final state fermion. Therefore, only the complete set of all QED diagrams 
is gauge-invariant. A 'nazve' subset of diagrams with ISR corrections corresponding to 
the diagrams with a (real or virtual) photon attached to the external electron legs is not 
gauge- invariant. This is different from the ZZ basic diagrams, where the electric charge 
flows continuously through the initial state. 

A straightforward solution of the problem would be a complete numerical calculation 
of all the 0{oi) corrections, including also intermediate and final state corrections. This 
was done for the on-shell case in |T^, but it seems to be incredibly complicated for the 
off-shell case (see |jl8|). In such a complete approach one must also include all electroweak 
corrections and properly treat the radiation from virtual intermediate VT^-states. 

We used a completely different approach to the problem, namely to take advantage of 
the fact that ISR corrections should yield the main fraction of the net correction. There- 
fore, we tried to define a gauge-invariant ISR correction by splitting the electrically neutral 
neutrino flow in the t-channel into two oppositely flowing charges -1 and +1. The charge 
-1 is then combined with the 'naive' ISR diagrams in order to build a continuous flow of 
electric charge in the initial state. The charge +1 is combined with the intermediate-state 
photon emission and is neglected here. This technique, which is explained in more detail 
in |T^, is called the Current Splitting Technique (CST). 

Within the CST, the ZZ and WW t-channel QED amplitudes are identical. The 
difference between the two cases arises at the level of cross-section calculations. For 
the WW case, there are st and tt interferences, while for the ZZ case one flnds t- and 
M-channel contributions and the tu interference. In all cases the cross-section has the 
following structure: 



j3 i3li?2,a / X 
" '^non-umvl'^/ 

dsids2ds' 



-pBi(Sl)pB2(s2) 



(27) 



with 



'5a(s, s'; Si, S2) 



l + S{s] 







+ (s ;si,S2), 



BiB2,a/ /. 



S ', Si, 82) 



(28) 



7Ya(s,s';si,S2) = H{s,s')aQ^^^'''{s';si,S2) 

+ af''^'^is,s';si,S2), (29) 

where a?^^^'^{s'; Si, S2) and a^^^'^'^{s, s'; Si, S2) are non-universal, non-factorizable soft 
and hard contributions. Equation (^) possesses several remarkable properties. The 
leading ISR correction contributions to the cross-section, containing mass singularities 
via Pe, factorize for any a = st, sm, tw-interferences and t- and -u-channel exchanges. This 
is necessary to ensure that the gauge cancellation is not spoiled. 
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The non-universal terms are calculated so far only for the lyiy-case |jT^. As one 
should expect, the non-universal contributions do not contain mass singularities. An 
analogous study for the ZZ case is in progress . 

The numerical influence of the universal part of the ISR corrections (i.e. setting non- 
universal parts equal to zero) is presented in figures 3-5 for the WW ^ ZZ^ and ZH cases. 
Universal ISR yields large, negative contributions in the vicinity of the threshold. At 
high energies these corrections are positive and one observes the effect of the radiative 
tail similar to the Z-peak. This radiative tail phenomenon is more pronounced than the 
high energy tail due to the bosons' off-shellness. We do not show the radiatively corrected 
on-shell cross-section 0"°^;^, but only mention that the relative differences between ag^j.^ 
and (Tg^j.^ on the one hand and between o"°JJjy and cr°^jy on the other are similar. 



Figure 7: The non-universal initial state and Coulomb corrections to off-shell VT-pair 
production. 

The effect of the non-universal contributions in the VTH^-case is illustrated in figure 7 
(solid line). This contribution is seen to be small, at LEP 2 energies it does not exceed 
0.4%. At high energies the relative contribution of the non-universal term becomes as 
large as 1.4% at = 1 TeV. To a great extent the smallness of the non-universal terms 
at high energies is due to the screening property of the non-universal corrections. They 
have a damping overall factor, 

st.t / / \ S\S2 /„„s 

a^'^(s;si,S2)~^. (30) 

The screening property ensures the unitary behaviour of the non-universal terms at high 
energy for the individual (st and t) contributions. 

In figure 7, we also show an important part of the final state corrections - the so- 
called Coulomb singularity. It yields a positive correction, which has its maximum value 
of about 6% at the threshold and vanishes at high energies. At 1/5= 1000 GeV it amounts 



to 0.75%. The Coulomb correction is taken into account according to equation (5) of pO |. 
However, at high energies, other final state corrections are important |1^ . 
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6 CONCLUSIONS AND PROSPECTS 



For a variety of purposes, the semi-analytical treatment of four-fermion production is an 
interesting alternative to the Monte Carlo approach. 
Within the GENTLE project, we have calculated so far 

(i) the 0{a) ISR corrections to and the average radiative energy loss (^^rad) in off-shell 
1^-pair production; 

(ii) the background contribution for this process with the simplest final state configu- 
ration; 

(iii) the off-shell ZH production with universal ISR corrections and the background 
contributions for the fi/ibb decay mode. 

We are presently studying 

(iv) the 0(a) ISR corrections to off-shell Z-pair production. 
We intend to study 

(v) background contributions to Z- and VF-pair production with other topologies; 

(vi) background contributions with t-channel exchanges and final states with and Ue 
or Pe- 

(vii) final state QED corrections to on-shell l^-pair production with the current-splitting 
technique; 

(viii) the inclusion of virtual weak corrections. 

The annihilation of two particles into four (five particles in the case of real 
bremsstrahlung corrections) has a limited variety of topologies in the tree approximation 
for the basic process. Thus, its systematic treatment is of principal theoretical interest. 
Apart from the case of LEP 2, it finds additional apphcations in e.g. the study of the Z 
line shape at LEP 1 (initial state QED or final state QED and QCD higher order pair 
production corrections) or in QCD corrections at the LHC. A study of certain problems 
beyond the Standard Model is, of course, also within reach. 
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